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Abstract 

Single-particle energies of A and S hyperons in several nuclei are ob- 
tained from the relevant self-energies. The latter are constructed within the 
framework of a perturbative many-body approach employing present realistic 
hyperon-nucleon interactions such as the models of the Jiilich and Nijmegen 
groups. The effects of the non-locality and energy-dependence of the self- 
energy on the bound states are investigated. It is also shown that, although 
the single-particle hyperon energies are well reproduced by local Woods-Saxon 
hyperon-nucleus potentials, the wave functions from the non-local self-energy 
are far more extended. Implications of this behavior on the mesonic weak 
decay of A hypernuclei are discussed. 

PACS Numbers: 21.80.+a, 13.75.Ev, 21.65.+f. 
Keywords: Hypernuclei, YN interaction, G-matrix, self-energy. 
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I. INTRODUCTION 



Hypernuclei are bound systems of neutrons, protons and one or more strange baryons, 
such as the A or S hyperons. Understanding the behavior of hypernuclei (how they are 
produced, their spectroscopy and decay mechanisms) has been the subject of intense inves- 
tigations during the last decades, see e.g., Refs. |I|-pi|. 

One of the main goals of such studies is to explore how the presence of the new degree 
of freedom (strangeness) alters and broadens the knowledge achieved from conventional 
nuclear physics. Several features of the A single-particle properties in the nucleus, being 
essentially different from those of the nucleon, have clearly emerged from these efforts. It 
is well accepted nowadays that the depth of the A-nucleus potential is around —30 MeV, 
which is 20 MeV less attractive than the corresponding nucleon-nucleus one. The spin-orbit 
splittings of single particle levels in A hypernuclei were found to be much smaller than their 
nucleonic counterparts typically more than one order of magnitude . Moreover, the A, 
contrary to the nucleon, maintains its single-particle character even for states well below 
the Fermi surface |T^JT3|| indicating a weaker interaction with other nucleons. Studies of 
the mesonic weak decay of light A hypernuclei [p!^-[T6[| have shown that the data |jT^ clearly 
favour A-nucleus potentials which show a repulsion at short distances. This seems also to be 
a characteristic of the S-nucleus potential for light S-hypernuclei [|T^ , which reproduces the 
recently measured bound S+ state in |;He with the in-flight '^B.e{K^ , tt^) reaction ||19|. This 
experiment confirms, with new and better statistics, the earlier results from the ^He(stopped 



K ,TT ) reaction 



Attempts to derive the hyperon properties in a nucleus have followed several approaches. 
A A-nucleus potential of Woods-Saxon type reproduces reasonably well the measured A 
single-particle energies of medium to heavy hypernuclei pT|-p3|. Non localities and den- 
sity dependent effects, included in non-relativistic Hartree-Fock calculations using Skyrme 
hyperon- nucleon (YN) interactions p^-p7|, improve the overall fit to the single-particle 
binding energies. The properties of hypernuclei have also been studied in a relativistic 
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framework, such as Dirac phenomenology p^j29[] or relativistic mean field theory p0|-|37[ . 



Microscopic calculations, which aim at relating the hypernuclear observables to the bare 
YN interaction, are also available. This approach is especially interesting because it can be 
used to put further constraints on the YN interactions, which are not completely determined 
by the limited amount of scattering data due, essentially, to the experimental difficulties 
associated with the short lifetime of hyperons and low intensity beam fiuxes. Microscopic 
hypernuclear structure calculations are performed with an effective YN interaction {G- 
matrix) obtained from the bare YN potential through a Bethe-Goldstone equation. The 
comparison with data may therefore further help in constraining the YN potentials. In 
several microscopic calculations a Gaussian parametrizations of the G-matrix calculated 
in nuclear matter at an average density [p7|-|40|| was employed. Furthermore, a G-matrix 
obtained directly in finite nuclei was used to study the single-particle energy levels of various 
hypernuclei or as an effective interaction in a calculation of the \''0 spectrum . 



Along similar lines, the authors of Ref. |^ derived microscopically the A self-energy in 
j^^O, starting from realistic hyperon-nucleon interactions. The starting point in the latter 
calculations is a nuclear matter G-matrix at a fixed energy and density, which is used to 
calculate the self-energy for the finite nucleus including corrections up to second order. The 
obtained self-energy for the A is non-local and depends on the energy of the hyperon. Solving 
the Schrodinger equation with this self-energy it is possible to determine the single-particle 
energies and wave functions of the bound hyperon. The approach also provides automati- 
cally the real and imaginary part of the hyperon optical potential at positive energies and, 
therefore, allows to study the hyperon-nucleus scattering properties. The method was first 



employed to study the nucleon and A properties in nuclei [^4|,^ and later applied to cal- 
culate the s-wave non-local A self-energy in ^^O, from which the single-particle A binding 
energy was obtained |^ . 



The aim of the present work is to extend the calculations of Ref. in order to derive 
s and p-wave A single-particle energies for a variety of A hypernuclei, from ^'^O to ^°^Pb. 
The S single-particle energies obtained for different potentials are also discussed. The YN 



potentials employed are the Nijmegen soft-core and the Jiilich interactions. Eventual 
differences in single-particle binding energies may therefore help in further constraining the 
YN interactions. 

In Sect. |T|, we present the formalism which consists of obtaining, first, the nuclear 
matter G-matrix at fixed density and starting energy parameter in the center of mass frame 
and, next, performing a second order calculation to derive the finite nucleus G-matrix from 
which one can obtain the hyperon self-energy for the different single-particle orbits. This 
self-energy can in turn be used to define a single-particle potential in a Schrodinger equation 
in order to obtain the corresponding single-particle binding energies. Nuclear matter results 



and a discussion on the convergence of our method is presented in Sect. [Ill A| . An important 
point raised in the latter section is that the single-particle hyperon energies calculated at 
first order with a nuclear matter G-matrix depend quite strongly on the fixed starting energy 
and density employed. The A single-particle energies for a variety of hypernuclei are shown 
in Sect, p.11 B| , where local equivalent Woods-Saxon potentials to represent our non-local 
self-energy are also derived. Finally, our conclusions are presented in Sect. |V]. 



II. FORMALISM 

In this section we present the formalism to obtain the hyperon single-particle energies 
in finite hypernuclei using an effective interaction derived microscopically from realistic 
YN interactions, for which we take the Nijmegen soft core [^] and Jiilich potentials. 
Although the formalism was already described in Ref. part of it will be repeated here 
in order to set up the notation used later in the description of the results. In the first 
place, we solve the G-matrix in nuclear matter at a fixed density, center of mass momentum 
and energy. Next, we transform this G- matrix into an effective interaction for the finite 
hypernucleus from which we can obtain the hyperon self-energy. Finally, we use this self- 
energy in a Schrodinger equation to derive the single-particle energies and the corresponding 
wave functions of the bound states. 
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A. Nuclear matter YN G-matrix 



The nuclear matter YN G-matrix is solved in momentum space and the two-particle 
YN states are defined in terms of relative and the center-of-mass momenta, k and K, given 

by 

Lr ^ 

Mn + My ' 
K = k^v + ky , 

where kjv and ky are the nucleon and hyperon momenta, respectively. Using an angle- 
averaged Pauli operator, we perform a partial wave decomposition of the Bethe-Goldstone 
equation, which, in terms of the quantum numbers of the relative and center-of-mass motion 
(RCM), is written as 

{{Y"N)kTKL{J)S'%\ G{u;nm) \{YN)klKL{J)ST,) = 

{{Y"N)k"l"KL{J)S'%\ V\{YN)klKL{J)ST,) 

+ E E E / k'^dk' {{Y"N)k"l"KL{J)S"T,\ V\{Y'N)k'l'KL{J)S%) 
I' S' y'=AS'^ 



g2 k'HMN+M ,) 



X {{Y'N)k'l'KL{J)S%\ G{u;nm) \{YN)klKL{J)ST,) , (1) 

where Q is the nuclear matter Pauli operator, V is the YN potential and ujnm is the nu- 
clear matter starting energy which corresponds to the sum of non-relativistic single-particle 
energies of the interacting nucleon and hyperon. Note that kinetic energies are used in the 
intermediate Y'N states and My — My accounts for the mass difference of the initial and 
intermediate hyperon. The variables A;, A;', k" and I, I', I" denote relative momenta and 
angular momenta, respectively, while K and L are the quantum numbers of the center-of- 
mass motion. The total angular momentum, spin and isospin projection of the YN pair are 
denoted hy J, S and T^, respectively. 
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B. The hyperon single-particle potential C/y in nuclear matter 



In the Brueckner-Hartree-Fock approach the hyperon single-particle potential Uy is ob- 
tained self-consistently by the following sum of diagonal G-matrix elements: 

UrikY) = I d'k^ {Yky, iVk^| G{eN{kN) + ey{ky)) {Yky, Nk^) , (2) 

where SM{Y){kN{Y)) = ^Ar{y)/(2MAr(y)) + UN{Y){kN{Y)) is the single-particle energy of the 
nucleon (hyperon). Using the partial wave decomposition of the G-matrix, the single-particle 
potential Uy can be rewritten as 

^y(M = ^ {2J + i){2T + i) 

z[zly -r i) J I 

X [''"^^^ k'^dkf{k,ky) {YN;klST,\G^{eN{kN) + ey{ky))\YN;klST,) , (3) 
Jo 

where an average over the hyperon spin and isospin (ty) has been performed and the weak 
center-of-mass dependence of the G-matrix has been neglected. In Eq. (^), A; is the relative 
momentum of the YN pair, ^y = M^/My, kmax is given by 

kp + ^yky 
i + 



and the weight function f{k, ky) by 



f{Kk 



Y) 



for k < '-^^ , 



for l^yky - {1 + ^y)k\ > kp , (5) 



In the discussion, we will associate the binding energy of the hyperon to its single-particle 
energy: 

By{ky)^e{ky) = ^ + Uy{ky) (6) 
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C. Evaluation of the hyperon self energy 



The self-energy of the A or E hyperon in a finite hypernucleus can be obtained, in the 
Hartree-Fock scheme, using a finite nucleus G-matrix as an effective YN interaction. How- 
ever, our G-matrix has been obtained in nuclear matter and, in particular, the intermediate 
propagator (that involves the Pauli operator and energy denominator) is very different than 
the corresponding one in a finite hypernucleus. One can, nevertheless, find the appropriate 
finite nucleus G- matrix, GpN) by relating it to the nuclear matter G-matrix through the 
following integral equation written in schematic form: 



GpN — G + G 
= G + G 
+G 



Q 



Q 



e ) FAT V e / NM 

9.) _ (9.) 

e J FN \ e J NM 

Q 



GpN 

G 



Q 

e / FN 



G 



9: 

e J FN 



e / NM 



G 



(7) 



e / NMi 

which involves the nuclear matter G-matrix (labelled G throughout the text) and the differ- 
ence between the finite nucleus and nuclear matter propagators. The latter account for the 
relevant intermediate states. In the practical calculations we will approximate the expansion 
up to the second order in the nuclear matter G-matrix 



GpN — G + G 



9: 

e J FN 



9 

e J NM. 



G 



(8) 



Therefore, in the evaluation of the hyperon self-energy we take into account the diagrams 
displayed in Fig. |I|, where the wiggly interaction lines represent the nuclear matter G-matrix. 
Diagram |I|(a) represents the first term on the right hand side of Eq. (P) which, by analogy 
to the nuclear case, will be refered to as the Hartree-Fock contribution. Diagram |l](b) stands 
for the second order correction, where the intermediate propagator has to be viewed as the 
difference of propagators appearing in Eq. (^). 

We will consider the incoming (outcoming) hyperon as a plane wave and the nucleon 
hole states as harmonic oscillator ones, so the two-body wave function will be a mixed 
representation of both single-particle states given by 



\{nhlhjhtz^){kYlY3YtzY)JTz) = J kldkhRnf^i^akh) \{khlhihtzh){kYlY3YtzY)JTz) , (9) 

where Uhlhjhtzh and kYWiYizy are the quantum numbers of the nucleon hole state and the 
hyperon state, respectively. Further, a is the oscillator parameter appropiate to describe 
the single-particle wave functions of the bound nucleons in the nuclear core. It is defined as 

a = , (10) 

with Tiuj chosen as the following function of the mass number 

huj = 45^-^ - 25A-i . (11) 

Typical matrix elements needed in the calculation are 

{{kYlY3YtzY){f^hlhihtzh)JTz\ G\{k'YlY3YtzY){'n'hlh3htzh)JTz) , (12) 
for the Hartree-Fock diagram of Fig. |l|(a), or 

{{kYlY3Ytzy){nhlh3htzJJT,\ G \{YN)klKL{J)ST,) , (13) 

appearing in the second-order diagram of Fig. |T|(b). 

Using Eq. (^, the mixed representation states appearing in Eqs. (|T^ and (|T^ can be 
expressed in terms of momentum and angular momentum variables in the laboratory frame. 
Then, with appropriate transformation coefficients [4S,|^, one can express the two-body 
states with laboratory coordinates in terms of the variables in the RCM system used in the 
solution of the G-matrixQ 



^Note the distinction between ka and k and la and With the notation ka or la we will refer to 
the quantum numbers of the single-particle state, whereas / or k without subscripts refer to the 
coordinates of the relative motion. 
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\{kalajatJ{hhjbt.,)JT,)= E I k'dk J K^dK 

ILXSJ 



la lb A 

lie 

2 2'-' 
ja 3b J 



L I A 



S J J 



(14) 



X {klKL\kalahlb) \klKL{J)SJT,) , 
where x = \/2x + 1 and {klKL\kalakblb) are the transformation coefficients from the RCM 
system to the laboratory system. 

We can construct now the expressions for the various diagrams considered in this work. 
The ffist order term of Fig. |l](a) yields a real and energy-independent contribution to the 
self-energy given by 

VHFikYkylYjYtzY) = T^E E 

X {{kYWjYtzY){nhkjhtzJ-JTz\ G\{kYlY3YtzY){'nhlh3htzh)JTz) , (15) 

where kY{kY)lYjYtzY are the quantum numbers of the incoming(outcoming) hyperon. 
The computation of the contribution coming from the two-particle-one- hole {2plh) diagram 
of Fig. |T](b) requires a little more work. First, we evaluate the imaginary part of the second 
term in Eq. (0). This term has an an explicit energy dependence. It reads 



y^tplhi^y^Y^YiYtzY^) 



jY niJhjht.,^ J ILSJY'=AT. 

X {{k'ylYjYtzY){nhlhJhtzJJTz\ G \{Y'N)klKL{J)SJTz) 
X {{Y'N)klKL{J)SJT,\ G\{kYlYjYtzY){nhlh3htz,)JTz) 



k^dk / K^dKJT 



X Tl5 [uo + Eh 



2{Mn + M, 



Y', 



2MnMy> 



My' + My 



(16) 



where uj is the energy of the hyperon measured with respect to the hyperon rest mass. The 
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single-hole energies Eh have been taken equal to the experimental single-particle energies in 
most of the nuclei studied (e.g., ^^C, ^^O, ^°Ca) and have been calculated from a Woods- 
Saxon potential with Spin-Orbit and Coulomb terms appropriately fitted in other cases (e.g., 
^°Zr and ^''^Pb). The quantities klKL{J)SJTz are the quantum numbers of the intermediate 
Y'N states. Next, we obtain the real part through a dispersion relation 

;(1) (U IJ 1 A ^ .,^ _ ^ [°° ^2plh{kYk'YlYjYtzY^') J, 



V'.Wvk'ylyJytzYiv) = - '^^P^nv^^'Y^rjr^ZY^; ^^,^ ^^^^ 

TT J-oo ijj' — UJ 

where P means a principal value integral. 

Finally, we must subtract the IjAh correction term coming from the nuclear matter 
intermediate propagator (third term in Eq. (H) ). It reads 

JY n^lhihU^ J lLSJY'=kT.-' 

X {{k'ylY3Ytzy){nhlh3ht-JJTz\ G \{Y' N)klKL{J)SJTz) 

X {{Y'N)klKL{J)SJT,\ G\{kylyJyt,^){nhkjhtz,)JTz) 

where Q is the nuclear matter Pauli operator and ujnm is the nuclear matter starting energy. 
This term only contributes to the real part of the hyperon self-energy and avoids the double 
counting over intermediate Y'N states contained already in the nuclear matter G-matrix of 
the Hartree-Fock contribution Vhf- 

In summary, the self-energy of the hyperon reads 

T,{kykylyjyu) = V {kykyly jyUj) + iW {kyk'ylyjyu) , (19) 

with the real part given by 

VikykylYjYUj) = VHF{kYk'ylYjY) + V2pih{kYk'ylY jyu) - Vf^i^ikykyly jy) (20) 

and the imaginary part by 

W{kyk'ylyjyu) = W^H^ikyk'ylyJyU). (21) 
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The self-energy can then be inserted as a single-particle potential in a Schrodinger equa- 
tion in order to investigate bound and scattering states of a hyperon in a finite nucleus. The 
different approximations to the self-energy, i.e., whether we include the 2plh contribution 
or not, result in different single-particle hamiltonians. We solve the Schrodinger equation 
by diagonalizing the corresponding single-particle hamiltonian in a complete basis within a 
spherical box of radius Rbox following the procedure outlined in This method is espe- 
cially suitable fon non-local potentials defined either in coordinate or in momentum space 

m. 



III. RESULTS AND DISCUSSION 



In this section we present and discuss results for A and S hypernuclei using two realistic 



interactions: Nijmegen soft core EH] and Jiilich B [B7 . 



A. Dependence of results on the starting energy 

As described before, our method provides the binding energies of the different hyperon 
orbits in finite hypernuclei starting from a G-matrix calculated in nuclear matter in the YN 
center-of-mass frame at fixed starting energy unm and Fermi momentum kp- By adding the 
2plh correction to the Hartree-Fock term one incorporates, up to second order in the nuclear 
matter G-matrix, the correct energy dependence and Pauli blocking factor in the finite 
nucleus. Therefore, the complete calculation {HF + 2plh) has to be viewed as a Hartree- 
Fock approach which uses a effective interaction derived microscopically with the appropriate 
density and energy dependence of the hypernucleus under study. This is in contrast to 
previous calculations p7|-^ where the determination of the finite hypernucleus effective 
interaction from the nuclear matter G-matrix implied a sort of average over the nuclear 
density. In these works several local and energy independent effective YN interactions of 
Gaussian form (YNG) were derived by parametrizing the corresponding nuclear matter G- 
matrices obtained from various YN potentials. The parametrization of the G-matrix into 
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a local effective interaction YNG to be used in finite hypernuclei calculations required the 
use of an appropiate value of the Fermi momentum kp. This value was determined, for each 
nucleus, by averaging the corresponding nuclear density weighted by the modulus squared of 
the A single-particle wave function of the single-particle level under study. The parameters 
of the effective YNG interaction were adjusted to reproduce the A potential energy f/A(0) 
in nuclear matter at the average value oikp- With these parametrizations, A single-particle 
energies and excited hypernuclear levels in several A hypernuclei were obtained through a 
shell-model calculation, with the aim of learning about the bare YN interaction. 

It therefore seems appropriate to explore, using our method, how much the hyperon 
single-particle energy depends on the starting energy (and density) of the nuclear matter G- 
matrix used in the calculation. This will allow us to assess how reliable energy independent 
effective interactions obtained from local density averages might be. 

Let us first show, in Fig. ^ the binding energy BY{k = 0) of a A (curves on the left) or 
a S (curves on the right) in nuclear matter at kp = 1.36 fm~^ as a function of the starting 
energy parameter u = ujnm + ^ =< > +BY{k = 0) + A, where < B^^ >= —50 MeV is 
an average of the nucleon binding energy over the Fermi sea and A = My — Ma. The long- 
dashed (full) lines are for the Nijmegen soft core (Jiilich B) interaction. An estimate of the 
self-consistent solution is obtained where the hne u = unm + A =< B^ > +BY{k = 0) + A 
crosses the calculated values of Byi^k = 0). This is indicated by the dotted hues in the 
figure. In the case of the Nijmegen interaction we obtain u = —74.3 MeV (-Ba(O) = —24.3 
MeV) for the A and u = 15.8 MeV (5s(0) = -11.7 MeV) for the S, whereas, in the case 
of the Jiilich interaction, u = -80.2 MeV for the A {Ba{0) = -30.2 MeV) and oo = -36.0 
MeV (Ss(0) = -63.5 MeV) for the E. 

Several features emerge from Fig. ^. First, the S hyperon is unrealistically overbound in 
nuclear matter by the Jiilich interaction. It is therefore necessary to readjust the parameters 
of this interaction if one wants to use it in shell model calculations of S hypernuclei. Secondly, 
we observe that the energy dependence of f/y(A; = 0) is slightly stronger in the case of the 
Nijmegen interaction, especially for the S hyperon which is more sensitive to the SA^ — AA^ 
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coupling because the starting energy is closer to the energies of the intermediate AA^ states 
(which propagate with the kinetic energy spectrum). Finally, we observe that the A binding 
energy varies at most by 10 MeV in a starting energy range of 80 MeV, while the variation 
of the S binding amounts to twice as much. As we will see below, this has consequences in 
the results for finite hypernuclei. 

In Tables | and || we show the binding energy of the A and respectively, in y O. The 
columns denoted by HF correspond to our lowest order calculation (see Eq. (|T5|) ) which 
uses, as effective interaction, the nuclear matter G-matrix calculated in the YN center of 
mass frame at fixed energy (shown in the first column) and density {kp = 1.36 fm~^). 
Columns labelled {HF + 2plh) include the 2plh corrections (see Eqs. (|1^) and (|18])) to 
bring the nuclear matter G-matrix to the finite nucleus one, with the proper energy and 
density dependence. We see that the lowest order results depend quite strongly on the 
starting energy used, especially in the case of the S hyperon as would be expected from the 
nuclear matter results shown in Fig. 0. However, it is worth noticing how, no matter what 
starting energy is used in solving the nuclear matter G-matrix, the corrected calculation 
[HF + 2plh) ends up giving practically the same result for the hyperon binding energy. 
Particularly stable are the results for the A hyperon. This weaker energy dependence, seen 
already in the nuclear matter results of Fig. H, is due to the fact that the energies involved 
in the calculation lie further away from the intermediate YN states, which propagate with 
kinetic energy, and therefore the strong AA^ — T,N coupling is less pronounced. 

In Tables |T| and |V| we show the binding energy of the A and S°, respectively, in y O 



using nuclear matter G-matrices calculated at several values of the Fermi momentum and 
a fixed value of the starting energy {u = —50 for the A and = for the S). The lowest 
order calculation for the hyperon single-particle energy, shown in the second column, depends 
quite strongly on the value oikp. However, one finds again that, when the 2plh correction is 
included to incorporate the proper intermediate propagator of the finite nucleus, the results 
nicely converge to practically the same value, no matter what was the density used in the 
solution of the nuclear matter G-matrix. 
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These results are interesting because they confirm that the finite nucleus Gi^Tv-niatrix 
is already well approximated by the second order term in the expansion in terms of the 
nuclear matter G-matrix, see the second and third terms in Eq. (|]). The correction, whose 
size depends on the starting energy or Fermi momentum used in the solution of the nuclear 
matter G-matrix, already leads to practically the same value for the hyperon single-particle 
energy. Higher order terms could only help in bringing the results closer than what they 
already are. Moreover, our results also show that in some cases the correction is quite 
appreciable, not only for the E° binding energies shown in Tables f| and |^ but also for the 



A energies in the case of the Nijmegen interaction. Therefore, if the correction is taken only 
approximately through an averaged nuclear matter G-matrix []37|-^ it may not lead to the 



proper effective interaction in the finite nucleus one is studying. These words of caution are 
particularly relevant in the case of the S hyperon where the corrections are very large. 

We note that the A single-particle energy obtained in the case of the Nijmegen soft 
core interaction is in excellent agreement with that obtained by Halderson (see column 3 
in Fig. 7 of Ref. |^T|), where the G-matrix was calculated directly in the finite nucleus 
for various Nijmegen interactions. Our method must be viewed as an alternative way of 
building up a finite nucleus effective interaction. It was already shown there that the Pauli 
corrections, which are a source of nucleus dependence, were very large for the Nijmegen soft 
core potential. This again supports our believe that calculations based on nuclear matter 
G-matrices at an average density will carry uncertainties tied to the chosen value of the 
Fermi momentum. In particular, the whole purpose of using these microscopic calculations 
to constrain the YN force cannot be achieved if a different value of /cj? is used for each YN 



interaction when studying the same hypernucleus ||3^ . 

In the studies of single-particle binding energies below, we will refrain from a study of 
the S binding energy since the results for ^^O give single-particle binding energies which are 
much too attractive. Several analysis of atomic data pO| , |5T| suggest a S well depth similar 
to that of the A p^EH] and more recent analysis did not discard S-nucleus potentials 



showing an inner repulsion. Moreover, E hypernuclear spectra from {K , vr ) reactions 
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suggest a relatively shallow S-nucleus potential Therefore, although the amount of data 
is limited, there is no experimental evidence for such strongly bound S hyperons. This 
clearly points to a weakness in the present YN interactions, hinting possibly at a too strong 
TjN — AN coupling in the interactions. 



B. A Single-particle states 

Once the method is well established and tested for the specific case of j^^O, it is the right 
moment to study the systematics of the A binding energy through the periodic table. To this 
end, values of the A single-particle binding energies obtained in what has been called HF and 
HF + 2plh approximations are reported in Table together with the available experimental 
data. These binding energies have been calculated using the energy-independent version with 
parameter set B of the Jiilich potential |^. The results for the Nijmegen soft core potential 
^ have not been considered in this section because the corresponding prediction for the A 



binding energy in nuclear matter is —23.4 MeV, about 7 MeV weaker than the prediction of 
the Jiilich model (—30.2 MeV) which agrees well with the extrapolated experimental values 
of the single-particle A binding energies with increasing mass number. 

The agreement with the experimental data is rather good. For convenience in the tech- 
nicalities of the algorithm we have always considered hypernuclei with a number of nucleons 
closing a subshell plus a A. Unfortunately, experimental data for those nuclei do not always 
exist and, as indicated in the Table 0, we have taken the closest representative nucleus 
for which the experimental information is available. Nevertheless, the differences between 
the calculated and the experimental values should not be associated to this fact but to the 
approximations used in the calculation or to the potential itself. 

For the density and starting energy used to calculate the A-nucleon G-matrix in nuclear 
matter, which has been used as effective interaction in our finite nucleus calculation, it turns 
out that the 2plh correction is always attractive. As discussed in the previous section, if we 
had used other starting values for the density or the starting energy we would have ended 
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up with different 2plh corrections but with the same value for HF + 2plh. 

In agreement with the experimental information, the difference between the P3/2 and pi/2 
A single-particle binding energies associated to these partial waves is very small. Note that 
the pi/2 energy is lower than the p3/2- This is a characteristic of the Jiilich interaction which 
yeilds too much attraction in the ^5*1 partial wave, as noted already in [§^,0 where the 
O"*" and 1"*" states of ^He were calculated and showed to appear in reverse order from the 
experimental values. 

The calculated A single-particle energies for ^°^Pb appear clearly overbound with respect 
to the experimental data. This is due to the fact that the distortion of the plane wave 
associated with the nucleon in the intermediate state of the 2plh diagram of Fig. |^(b), 
necessary to ensure its orthogonalization to the nucleon hole states, has been taken only 



approximately. The orthogonalization procedure is described in Ref. [44] and has been 



optimized for the case of \ 0. Actually, this feature is already sizable for ^ Zr and in the 
case of ^°^Pb leads to a result which is more bound than a A in nuclear matter. 

Traditionally, the systematics of A single-particle binding energies has been studied by 
using a phenomenological Woods-Saxon potential 

^^^(^) = TT m/v (22) 

1 -I- exp [{r — R)/a\ 

with a fixed diffusivity a and depth Vq, and a radius R = tqA^^^. A good parameterization 
of the experimental data is obtained with Vq = —30.7 MeV, tq = 1.1 fm and a = 0.6 fm 
25| . A more refined analysis of the same authors allows for a smooth dependence of tq 



in A and a slightly shallower potential (Vq = —28 MeV) with a larger radius, ro(A) = 
(1.128 + 0.439yl~^/^) fm, provides a better agreement with the experimental A binding 
energies. We have performed a similar analysis for the calculated A binding energies. In 
principle, the calculated self-energy is non-local both in fc-space and in r-space. However, 
in a previous work |^ we have shown that one can generate a local representation of 
the self-energy by performing an appropriate average of the non-local self-energy EQ,(r, r'), 
where a indicates the quantum numbers of the single-particle state, over the coordinate r'. 
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This local representation might, in first approximation, be characterized by the shape of 
a Woods-Saxon potential. Instead of doing this average, a possible strategy is to assume 
a Woods-Saxon shape, fix the depth and the diffusivity independent of the mass number 
and determine the radius R by requiring the Woods-Saxon potential to reproduce the same 
eigenvalue than the microscopic non-local energy dependent self-energy. A reasonable value 
for the depth Vq is the A binding energy in nuclear matter, which is taken to be —30.2 MeV, 
and for the diffusivity a = 0.6 fm. The resulting values of R when we apply this procedure 
to the s deepest state of j^C, and ^^Ca are 2.25 fm, 2.53 fm and 3.82 fm respectively. 
Fitting these three values with a functional form similar to the one used in Ref. for the 
analysis of the experimental data one obtains ro{A) = (1.229 - 1.390^-^/3) fm. In order to 
visualize the quality of these fits, we show in Fig. ^ the binding energies for the s and p waves 
of j^C, ]^0, j^Ca. and ^^Zr calculated with our non-local self-energies (triangles) together 
with the values obtained with a Woods-Saxon potential with the parameters just defined 
above (solid lines). As the spin-orbit splitting is so small we have reported the average 
value of the P3/2 and pi/2 energies obtained from the non-local self-energies and have not 
considered any spin-orbit term in the adjusted Woods-Saxon potential. The results of j^^Ph 
have not been included in the plot because, as mentioned before, the s wave binding energy 
was larger than the binding energy in nuclear matter which we have taken as the depth 
of the Woods-Saxon potential. The calculated binding energies are well reproduced by the 
Woods-Saxon shape and, as expected, both partial waves extrapolate to the binding energy 
for nuclear matter. 

Of course the binding energies are not enough to characterize the single-particle states 
since potentials giving rise to the same binding energies can generate substantial differences 
in the corresponding wave functions. Therefore, in order to analyze the microscopically 
calculated self-energy we should also study the single-particle wave functions. 

To have a measure of the goodness of the wave functions generated by the Woods- 
Saxon potential, we calculate their overlap with the wave functions obtained by solving the 
Schrodinger equation using the self-energy. The overlaps for ]^C, ^O, j^Ca and ^Zi are 
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0.9917, 0.9869, 0.9924 and 0.9853 respectively, which are not close enough to 1 to guarantee 
the equality of the wave functions. This is visualized in Fig. ^ , where the wave function for 
the s wave in j^^O obtained with the Woods-Saxon potential (dashed-line) is compared with 
the one obtained directly from the self-energy (solid line). 

Another possibility would be to keep the diffusivity fixed and adjust the radius R and the 
depth Vq to reproduce the eigenvalue and to maximize the overlap with the eigenfunction 
provided by the self-energy. The values of Vq by applying this procedure to the s wave of 
the three lighter nuclei considered above are respectively —23.11, —23.56 and —27.84 MeV 
whereas the values of the radius R are 2.92, 3.32 and 4.39 fm. With these values of Vq and R 
the overlaps are 0.9999 for the three nuclei considered. The eigenfunction obtained by this 
procedure for the is also drawn in Fig. ^ (dot-dashed line) and shows a large overlap 
with the self-energy eigenfunction (solid-line). 

In conclusion, the single-particle energies of closed-shell nuclei with one A are well repro- 
duced by using both the microscopic self-energy or the simpler parametrization of Woods- 
Saxon type in the Schrodinger equation. However, the wave functions provided by the micro- 
scopic self-energy differ from the ones originated by a Woods-Saxon with a fixed depth and 
diffusivity and a A-dependent radius. It is important to note that the mean square radius of 
the self-energy eigenfunction is larger than that from the corresponding Woods-Saxon wave 
function. This can have important consequences in the study the mesonic decay of these A 
hypernuclei. Indeed, it has been observed that the mesonic decay rates of light hypernuclei, 
such as ^H, ^He and \iie, could be better reproduced if the A wave function was pushed 
out to the surface by the effect of a repulsive hyperon-nucleus potential at short distances. 
This would favour the mesonic decay of these hypernuclei because the A would be exploring 
smaller nuclear density regions and the Pauli blocking effects, which prevent the mesonic 
decay from occuring, would be less pronounced. The mesonic decay rates of light hypernuclei 
have been calculated using repulsive A-nucleus potentials at short distances that have been 
obtained either phenomenologically []TB[, from a quark based bare YN interaction []T3[ or 



from a microscopic YNG effective interaction folded with an extremely compact He density 



1^. At present, no calculation exists that combines the use of an YN effective interaction 



with an appropiate density treatment of the host nucleus. Our method provides such ingre- 
dients and has been shown to produce A wave functions that are pushed out to the surface. 
This is a consequence of the non-localities of the self-energy and might not be related to a 
repulsive character of the A-nucleus potential at short distances. The implications of our 
results on the mesonic decay of A hypernuclei will be explored in a future work. 



IV. CONCLUSIONS 

We have analyzed a method to obtain the effective hyperon-nucleon interaction in finite 
nuclei based on an expansion over a G-matrix calculated in nuclear matter at fixed density 
and starting energy. The purpose of this study is to set up a reliable frame for hypernuclear 
structure calculations with the aim of obtaining information about the hyperon-nucleon 
interaction, complementary to that provided by hyperon-nucleon scattering experiments. 

We have shown, by explicit calculation of the A and S single-particle energies in yO, 
that truncating the expansion over the nuclear matter G-matrix at second order gives re- 
sults that are very stable against variations of the density and starting energy used in the 
G-matrix. Moreover, both first and second order terms depend quite strongly on those pa- 
rameters. This is an indication that the density dependent effects considered when treating 
explicitly the finite size of the nucleus are very important and, therefore, they might not be 
well approximated by energy independent and local effective interactions which start from 
a parametrized nuclear matter G-matrices evaluated at an average density. We note that 
the use of local and density averaged effective interactions can be extremely useful in de- 
tecting similarities and differences among the various hyperon-nucleon potentials. However, 
if the aim is to fine-tune the bare YN interactions to reproduce the spectroscopic data of 
hypernuclei, an appropiate effective interaction for the hypernucleus under study, as the one 
provided by our method, is in order. In particular, the two interactions employed in the 
present work give rise to very attractive S binding energies, while the S~ atomic data and 
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{K~,7r~) spectra seem to indicate E-nucleus potential depths similar to that of the A or 
even shallower. 

Although the method can be viewed as an alternative way of building up a finite nucleus 
effective interaction, it provides also the complete energy dependence of the hyperon self- 
energy. This allows in turn for a study of not only the bound states, as done here, but also 
the scattering states. This is especially of interest in the analysis of hypernuclear production 
reactions which yield a large amount of quasifree hyperons. 

We have obtained local Woods-Saxon A-nucleus potentials that reproduce the A single- 
particle energies of several hypernuclei. However, the wave functions obtained from our 
non-local self-energy are far more extended and can be simulated only when we allow the 
Woods-Saxon potential to have an 74-dependent depth and a relatively larger radius. This 
can have important implications on the weak mesonic decay of A hypernuclei, which will be 
explored in a future work. 
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TABLES 



TABLE L Dependence of the A single-particle energy in Jy O on the starting energy of the 
nuclear matter G-matrix. Our notation is lo =< Bn > +B\{k = 0), with < Bn >= —50 MeV. 



kp = 1.36 fm-i 


Nijmegen 


Jiilich 


UJ 


HF 


HF + 2plh 


HF 


HF + 2plh 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


-100 


-3.83 


-7.43 


-9.25 


-11.85 


-80 


-4.76 


-7.39 


-10.15 


-11.83 


-50 


-5.59 


-7.36 


-11.73 


-11.84 



TABLE IL Dependence of the TP single-particle energy in ^oO on the starting energy of the 
nuclear matter G-matrix. Our notation is oj =< Bjq > +Bj:{k = 0) + A, with < Bn >= —50 
MeV and A = Ms - Ma. 



kp = 1.36 fm-i 


Nijmegen 


Jiilich 


UJ 


HF 


HF + 2plh 


HF 


HF + 2plh 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 





-0.16 


-22.79 


-36.70 


-50.70 


20 


-2.01 


-23.35 


-40.38 


-50.94 


50 


-10.65 


-24.62 


-51.34 


-50.38 
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TABLE III. Dependence of the A single-particle energy in y^^O on the Fermi momentum of the 
nuclear matter G-matrix. Our notation is u =< Bn > +B\{k = 0), with < Bn >= —50 MeV. 



u = -50 MeV 


Nijmegen 


Jiilich 


kp 


HF 


HF + 2plh 


HF 


HF + 2plh 


(fm-i) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


1.00 


-9.33 


-7.30 


-13.71 


-11.74 


1.25 


-7.66 


-7.34 


-12.56 


-11.83 


1.36 


-5.59 


-7.36 


-11.73 


-11.84 



TABLE IV. Dependence of the S single-particle energy in ^oO on the Fermi momentum of 
the nuclear matter G-matrix. Our notation is uj =< Bn > +Bj^{k = 0) -|- A, with < Bn >= —50 
MeV and A = Ms - Ma. 



a; = MeV 


Nijmegen 


Jiilich 


kp 


HF 


HF + 2plh 


HF 


HF + 2plh 


(fm-i) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


1.00 


-2.64 


-25.03 


-43.94 


-51.74 


1.25 


-0.82 


-23.37 


-39.41 


-51.16 


1.36 


-0.16 


-22.79 


-36.70 


-50.70 
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TABLE V. A binding energies in the lsi/2) lP3/2 ^^id lpi/2 single-particle orbits for different 
nuclei. The available experimental data, indicating the hypernucleus for which they have been 



measured, are taken from the compilation of M supplemented by new measures reported in |12] 



and |13|. All the results have been derived from the Jiilich B interaction. 



Hypernuclei 


Orbit 


HF 


HF + 2plh 


Exp 


13p 


l'Sl/2 


-7.93 


-9.48 


-11.69 


A ^ 


l'Sl/2 


-10.15 


-11.83 


(a 0) 
-12.5 




1^3/2 
lPl/2 


-0.08 


-0.87 
-1.06 


-2.5 (Ip) 


A Ca 


lSl/2 


-16.85 


-19.60 


(a Ca) 
-20. 




lP3/2 
lPl/2 


— D. ( u 

-6.92 


-9.92 


— IZ. {Lp) 




l'Sl/2 
lP3/2 
lPl/2 


-22.24 
-14.74 
-14.86 


-25.80 
-18.19 
-18.30 


(fZr) 
-23. 
-16. (lp) 




1^1/2 


-26.28 


-31.36 


(rpb) 

-27. 




l]33/2 
lPl/2 


-21.22 
-21.30 


-27.13 
-27.18 


-22. (\p) 
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FIGURES 




(a) (b) 

FIG. 1. Diagrams through second order in the interaction YN G (wavy line) included in the 
evaluation of the hyperon self-energy. Diagram (a) is the first order term, while (b) is the second 
order 2plh correction. 
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FIG. 2. Dependence of the hyperon binding energy Byik = 0) in nuclear matter on the starting 
energy lu. The curves on the left are for the A, whereas the ones on the right are for the S. 
Long-dashed (full) lines correspond to Nijmegen soft core (Jiilich B) interaction. The dotted lines 
show the position of the self-consistent solution for Byik = 0). 
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FIG. 3. Calculated A binding energies in Is and Ip single-partkce orbits for different nuclei. 
The curves correspond to the solutions obtained for a Woods-Saxon potential whose parameters 
are defined in the text. 
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FIG. 4. The wave function in r-space for the lsi/2 A in j!JO obtained from the A self-energy 
(soUd hne) is compared with the ones obtained using a Woods-Saxon potential of fixed 
(j4-independent) depth (dashed line) or with both radius and depth adjusted (dot-dashed line) 
to maximize the overlap with the wave function provided by the self-energy. 
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